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Abstract 



In this letter, we present the general form of equations that gener- 
^ ' ate a volume-preserving flow on a symplectic manifold {A4, lo). It is 

shown that every volume-preserving flow has some 2-forms acting the 
role of the Hamiltonian functions in the Hamiltonian mechanics and 
the ordinary Hamilton equations are included as a special case with a 
2-form ^^^^ H 00 where H is the corresponding Hamiltonian. 

In classical mechanics, most of the systems can be described by the ordi- 
nary Hamilton equations 

^ dpi ' dq'^ ' 
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For an autonomous system, namely, a system with the Hamiltonian man- 
ifestly independent of the time t, the Hamiltonian if is a function on the 
phase space Ai which is precisely the cotangent bundle of the configuration 
space. Abstractly, a Hamiltonian system can be defined on a symplectic man- 
ifold {Ai,uj) where is a 2n-dimensional orientable different iable manifold 
endowed with a symplectic 2-form uj at every point of Ai. For the details 
of Hamiltonian mechanics on a symplectic manifold, see, for example, 
Then, according to Darboux's theorem, for each point x of A^, there is a 
coordinate neighborhood {U; q^, . . . , . . . ,Pn) of x such thatQ 

uj\u = dpi A dq\ (2) 

Thus a Hamiltonian system with the Hamiltonian H can still be described 
locally by eqs. (|l|). The solutions of eqs. (P are the integral curves of the 
Hamiltonian vector field 

■= —— - —— (3) 

dpi dq^ dq^ dpi 

on A4. It is well known that the above ordinary Hamiltonian system pre- 
serves both the symplectic 2-form u and the volume form ^ cj", where 
n = I dim and is the fc-fold wedge product of uj. 

It is well known that the ordinary Hamilton equations can be expressed 

as 

- zxh^ = dH, (4) 

where i denotes the contraction. It is easy to show the symplectic structure- 
preserving and volume-preserving properties for these ordinary Hamiltonian 
systems. Namely, 

jCx„io = (5) 

and 

/:x,K) = n(£x,^)Aa;"-i = 0, (6) 

respectively. 

However, there exist certain important mechanical systems that preserve 
the volume of the phase space only and cannot be transformed into the 
ordinary Hamiltonian systems. For example, the linear systems on 

q ' = -h, (7) 

^ Repeated indices are always assumed to be summed over 1 to n unless otherwise 
stated. For convenience, symbols such like are always abbreviated as 
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with constant coefficients kij are such kind of mechanical systems. Obviously, 
eqs. (1^ can be turned into the form: 

1 P« = -^-%^^ (8) 

with 

H=U'^p,p, + U,,q'q\ (9) 



I , . 1 



Sij 



- {kij + fcjj), aij — - {kij kji). (10) 



When one of aij is nonzero^, the system (^) is not an ordinary Hamiltonian 
system on M^". In fact, the corresponding vector field 

. d 

X = Xh - aij q^ — 

dpi 

of (^ is not even a symplectic vector field since the Lie derivative 

Cx.oJ = aij dg* A dq^ 

does not vanish in this case. On the other hand, the linear system (H) always 
preserves the volume form of the phase space because £x('^") = 0, as what 
can be easily verified. 

In this letter, we present a general form of equations that generate a 
volume-preserving flow on a symplectic manifold {A4,lj). It is shown that 
every volume-preserving flow has some 2-forms acting the role of the Hamil- 
tonian functions in the Hamiltonian mechanics and the ordinary Hamiltonian 
equations are included as a special case with a 2-form ^j^^- H uo where H is 
the corresponding Hamiltonian. 

^ Here is an example showing that even a conservative system can be transformed to 
be a non-Hamiltonian system. Consider a system consisting of two 1-dimensional Unearly 
coupled oscillators: 

miq^ ^ ~k{q^ - q^), m2 q"^ ^ {q^ - q^) . 



Obviously, such a system satisfies Newton's laws. Let kn = —k/mi, fci2 = k/mi, k2i = 
k/m2 and k22 = —k/m2. Then it is a system as described by eqs. (0), with the indices 
running from 1 to 2. When mi ^ m2, we have a system satisfying ki2 ^ fei- 
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The main idea for presenting such general equations for volume-preserving 
systems is based on the linear isomorphism between the space X^{Ai,uj) 
of certain volume-preserving vector fields and the space B'^"'~^{Ai) of exact 
(2n — l)-forms on Ai. This is, in fact, a direct generalization for the ordinary 
Hamiltonian system {n — 1). Namely, the linear isomorphism between the 
space X^^^{A4,u) of certain symplectic structure- preserving vector fields 
and the space B^{Ai) of exact 1-forms on A4. 

Briefiy speaking, a volume-preserving system, described by a volume- 
preserving vector field with respect to which the Lie derivative of a;" being 
zero, corresponds to a unique closed (2n — l)-form on M. and vise versa. 
Such a system that corresponds to an exact {2n — l)-form can be determined 
then by a 2-form (may up to a closed 2-form) , as if a Hamiltonian system can 
be determined by a function (the Hamiltonian) up to a constant function. 

In this letter, the main results will be quoted and outlined without de- 
tailed proof. For the details and extended topics, we refer to our forthcoming 
papers. 

For each integer k = 1, . . . ,n where n = ^ dimTW ^ 1, we may define the 
set of the symplectic vector fields and that of the Hamiltonian vector fields 
as 

X^{M,u;) = {Xe X{M)\C^{u^) =0}, 
X^{M,(jj) = {X e A'(A<) I -ix(a;^) is exact}, 

respectively. Here and hereafter X{A4) denotes the space of smooth vector 
fields on M.. Then it is easy to prove that 

Xi{M,uj) C ... C Xi{M,Lu) C X^+\M,uj) C ... C X^{M,uj), 
X^{M,uj) C . . . C X^{M,uj) C X^+\M,uj) C . . . C X^{M,uj) 

and 

In fact, Xg{Ai,uj) is a Lie algebra under the commutation bracket of vector 
fields, and X-^{A4,uj) is an ideal of X^{A4,u) because 

[Xi(M,u;),Xi(M,u;)] C X^{M,u;). 
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Let us consider the quotient Lie algebra 

H^^{M,uj) := Xi{M,u;)/X^{M,u;). (11) 

It can be called the fc-th Euler- Lagrange cohomology group. It is called a 
group because its Lie algebra structure is trivial: For each k = 1, . . ., n, 
H^i^{J^,uj) is an abelien Lie algebra. For n = 1, H^i^{J^,u) is the Euler- 
Lagrange cohomology first introduced in |^. And it can be proved that 
H^i^{Ai,uj) is linearly isomorphic to the first de Rham cohomology group 
HUM). 

For k = n, X^{A4,uj) is the Lie algebra of volume-preserving vector 
fields, including all the other Lie algebras mentioned in the above as its 
Lie subalgebras. Similarly to the case of k = 1, the n-th Euler-Lagrange 
cohomology group H^^{A4,uj) is isomorphic to the (2n — l)-th de Rham 
cohomology group H'^^^^ {A4) . To prove this, let us first introduce a lemma: 

Lemma 1 For each integer A; = 1, . . . , n (n ^ 1) and x G A4, the linear 
map Uk : T^M. — > K2k-i{T*JVi) defined by i^fc(X) = —ixi^^'') is injective, 
namely, ix(t^'^) = iffX. = 0. 

Since dimT^TW = dim A2„_i(T^A^), lemma |1| implies that z/„ : T^J^ — *• 
^2n--i{Tx-M) is a linear isomorphism, which induces a linear isomorphism 
i^n:X{M) — > n^''~\M) defined by 

z/„(X) = -ZxK) = -n (zx^) A u;--\ (12) 

Hence we obtain 

Theorem 1 The linear map Un : — > Q'^^''^^{Ai) is an isomor- 
phism. Under this isomorphism, X^{Ai,uj) and X{l{Ai,uj) are isomorphic 
to Z'^"'~^{M) and B'^"^^ (M) , respectively. 

In the above theorem, Z'^^~^{Ai) is the space of closed {2n — l)-forms on Ai 
and i?^""^(A4) is the space of exact (2n — l)-forms on M.. As a consequence, 
we obtain the following corollary: 

Corollary 1 The n-th Euler-Lagrange cohomology group H^^{Ai,uj) is lin- 
early isomorphic to H^^^ (Ai) , the {2n — l)-th de Rham cohomology group. 
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When is closed (i.e., compact and without boundary), such as a closed 
Riemann surface, H^i^{A4,uj) is linearly isomorphic to the dual space of 
H^i^iM), because H^^{M) = (Hl^-^M))* for such a manifold (see, for 
example, 0]). If is not compact, this relation cannot be assured. 

When n > 2, there is another lemma: 

Lemma 2 If n > 2, then, for an arbitrary integer k = 1, . . . ,n — 2 and a 
point X e M, a e A2iT*M) satisfies a Au'' = iff a = 0. 

Again this induces injective linear maps 

a I — > a Au^ 
for = 1, . . . , 72 — 2. Thus we can obtain that 

When k = n — 2, it can be easily verified from lemma that 

a ^ «Acu"-2 (13) 

is a linear isomorphism. If ri = 2, we use the convention that cu*^ = 1, namely, 
i = id : a i — > a. Then we can define a linear map making the following 
diagram commutative: 



d 



(14) 



Equivalently, given a 2-form 

" = ^ Q^3 dg* A dq' + A) dpi A dq' + \ P'^ dpi A dp,- (15) 

where Qij and P'-^ satisfy 

Q,i = -Q,„ P'' = -P'\ (16) 
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the vector field 

0(q;) = ^ o d o = u^^{da A a;""^) 

is in X^{M,uj). 

For convenience, we set 

X = n(n - 1) 0(q;) = n{n - 1) i/~^(dQ; A a;"-^), 

namely, ixl'^") = —n{n — 1) (da) Aa;"~^. It is easy to obtain that 

\ dq^ dpi dpj J dq^ \ dpj dq^ dq^ J dpi ' 



In fact, 

(ii[a) = da A a; 



.n-2 



,n-2 



1 do ■ ■ 1 P*-* 

+ ^ + ^) dp, A d?^- A dq'^ A c"-^ 

By virtue of the following two equations 

dpi A dp,- A dg*^ A o;"-^ = ^^dp^Au;"-^ ^ dp. A (18) 

n — 1 n — 1 

dpi A dg^' A dg'^ A u;'*-^ = dg'^ A u;""^ - dg^' A (19) 

n — 1 n — 1 

we can write dt(Q;) as 



n — 1 V 9g* ^g-? 9pj 
7 



Aa;'^-^ 



n — 1 V Sg-J dpj * 



1 /Mj 9 A] dQ 



2 e 



n(n — 1) V dq^ dpj J 



n(n — 1) V dq^ dpi dpj J 8q 



Comparing it with 



d^(«) = -T^, ^nO^) = ^xK), (20) 

n[n — 1) n[n — 1) 



we can obtain the expression of X, as shown in eq. 

Since both l and z/„ are hnear isomorphisms, we can see from the com- 
mutative diagram (|l^ that for each 2-form a on as in eq. (^5|) the vector 
field X in eq. (0) belongs to X{^{M,uj) and for each X G X{^{M,uj) there 
exists the 2-form a on satisfying eqs. But there may be several 

2-forms that are mapped to the same vector field X. For example, the vector 
field X in ([1^) is invariant under the transformation 

a I — >a + e (21) 

where 6' is a closed 2-form. 

Note that for the vector field X G X^{M., cj), the 2-form a is a globally 
defined differential form on A^. If X G X§{M.^uj)^ such a 2-form cannot 
be found if H'^^{A4.) is nontrivial. In this case, a can be still found as 
a locally defined 2-form, according to Poincare's lemma. Then the relation 
between X and the locally defined 2-form a, eq. is valid only on a 

certain open subset of A^, and, on the intersection of two such open subsets, 
the corresponding 2-forms are not identical. That is, the transformation 
relation of Qij, A^j as well as P*-' is not that of tensors on Ai. Instead, a 
transformation such like (pT]) or more complicated should be applied. 

No matter whether X belongs to X{^{Ai,uj) or X^{Ai,uj), namely, whether 
the 2-form a is globally or locally defined, the flow of X can be always ob- 
tained provided that the general solution of the following equations can be 
solved: 

_ dP'^ dAi _dAi 
^ dq^ dpi dpj ' 
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dpj dq^ dq^ 



= ^-^ + lp7- (22) 



This is just the general form of the equations of a volume-preserving mechan- 
ical system on a symplectic manifold {Ai,uj). 
For a function H on A4, let the 2-form be 

a = ^—Huj, (23) 
n — 1 

then eqs. (0) turn to be the ordinary Hamilton equations (p. Thus the 
ordinary Hamilton equations have been included as a special what 
is expected. 

Now suppose the symplectic manifold is M^" with the standard symplectic 
form u = dpi A dq\ When Qij = —aijpkq'', A} = H Sj and P*-' = with 
the constants aij and the function H as shown in eqs. (^ and ([lO|) , namely, 

a = ^^-^ Hu-^pkq'' aij dq' A dg^' , (24) 

then the equations (^) will turn into eqs. (^). Therefore the 2-form a in 
eq. ( p4D is one of general forms corresponding to a linear system, but it is 
not the unique. 

It is interesting that when if = in the above equation, all the coordi- 
nates are the first integrals of the linear system. Then all the canonical 
momenta pi = aij qH + pio where pio are constants. This can generalize to an 
arbitrary symplectic manifold for the equations (^), even though the sys- 
tem is no longer a linear system: If, on a Darboux coordinate neighborhood 
{U; q,p) in a symplectic manifold {J^,uj), the 2-form a satisfies 



lag 

9Pi 



A) dq^ + P^^ dpj = (25) 



for each i = 1, . . . ,n, then A'j = and P*-' = 0. Hence, on that coordinate 
neighborhood U, all the are constant. 
If we define a function tr a as 

aAa;"-i = — cu" (26) 
n 

for each 2-form a, then using the formula 

dpi A dq^ A cj"-^ = ^ cj", 
n 



9 



we obtain that 

tra = A^. (27) 

The above expression is obviously independent of the choice of the Darboux 
coordinates. Let Xtra be the Hamiltonian vector field corresponding to the 
function tra. Eq. (^) indicates that 

X = Xt,, + X' (28) 

where X' is the extra part on the right hand side of eq. (0), corresponding 
to the 2-form 

tr a , , 

a uj. 29 

n-l ^ ' 

If f{q,p) is a function on Ai, then the derivative / = ^f{qit),p{t)) along 
any one of the integral curves of eqs. ([2^ ) satisfies the equation 



/cu" = n(n- l)da Ad/ Acu"-2. (30) 
In fact, fit) = {C:}^f){q(t),p(t)). And, since X is volume-preserving, 
(/:x/)^" = C^ifio^) = dt^{fu-) + tj,d{fu-) = d(/^xc^"). 
Then according to eq. (|20|) , 



(£x/)cu" = -n(n-l)d(/di(a)) = -n(r2-l)d/Ad(aAcu"-2) 

= -n{n - 1) d/ A da A cu""^ = n{n - 1) da A d/ A cu"-^. 

Thus eq. (^) has been proved. 

Especially, when a = to, the system (^2]) turns out to be the usual 
Hamiltonian system, as we have mentioned. For such a Hamiltonian system, 
on the one hand, we can use eq. (|30| ) to obtain 

/ = n d{Huj) A d/ A a;"-2 = ndH Adf A a;""^ = tr {dH A df) w", 

namely, 

f = tT{dHAdf). (31) 
On the other hand, / can be expressed in terms of the Poisson bracket: 

df OH df dH 



f = {f,H} :=Xh/ 



(9g* dpi dpi (9g* 
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So we obtain the relation between the Poisson bracket and the trace of 2- 
forms: 

{f,H} = -tT{df AdH). (32) 

In order to develop a volume-preserving algorithm, Feng and Shang in 
presented a lemma as follows: The volume-preserving vector field X = 
{X^, . . . ,X"')'^ on can always be expressed by an antisymmetric tensor 
a'^ on aspl 

- It- (»' 

where are the standard coordinates on R". In fact, this lemma is not 
difficult to understand: A vector field on M", known as a Euclidean space or a 
Riemannian manifold, is volume-preserving if and only if it is divergence-free. 
Using the Hodge theory, the divergence of the vector field X reads divX = 
—SX = *d* X, where 6 is the codifferential operator and X = 6ij X^ dx^ is 
the 1-form, or the covariant vector field, obtained from X by using Einstein's 
method of lowering indices. Hence that X is volume-preserving, namely, 
divX = *d * X = 0, implies *X is a closed (n — l)-form. Then, according 
to Poincare's lemma, there is a 2-form a, say, such that *X = d* a. If a is 
assumed to be | aij dx* Ada;-' with aij = —aji = SikSji a^^ = a*-', we can obtain 
that X = (— l)""-*^ * d * a = 5ij ^^dx\ Then the vector field X satisfies 



eq. m- 

It is worthy to see that the formula ( P5| ) is quite similar to eq. (JT^) , only 
with some slight differences: 1. This formula is a statement on the Euclidean 
space with an arbitrary dimension n while eq. (|T7p is for a symplectic 
manifold, of 2n-dimensional. 2. It can be generalized to an n-dimensional 
Riemannian or pseudo- Riemannian manifold provided that the (n — l)-th 
de Rham cohomology group is trivial, with the ordinary derivatives being 
replaced by the covariant derivatives. As for eq. (0), when H^^~^{A4) ^ 
0, not every volume-preserving vector field satisfies it. This is similar to 
eq. (0). But a covariant derivative is not necessary in eq. (|1^. 3. When 
the symplectic manifold M. = M^", we can set 

, , _ / l\fQ-A^\fO -I\_fp-A 

la,,j2nx2n "(j P [l ~ [ A'^ Q 



^ In the quotation of this lemma, we have adopted the notations so as to accommodate 
ours in this letter. 
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where Q = {Qij)'l=,^j=„ P = (P'^)r=i,j=i, A = (4)r=i,,=i and J is the n x n 
unit matrix. Accordingly, we set {x^, . . . , x", x""*""^, . . . , x^") = {q^, . . . , q"',Pi, ■ 
Then eq. (^) reads 

dP'^ OA] d4 dQ,j 



^ dq^ dpj ' dq^ ~'~ dpj ' ' ^^^^ 



with the index i running form 1 to n. Comparing it with eq. ([17|), it seems 
that a condition 

tra = (35) 

could be imposed on the 2-form a in eq. {\L7\). 

Finally, it should be mentioned that the volume-preserving systems play 
an important role in physics, especially in statistical physics. However, this 
has not been noticed widely yet. In this letter we have not explained why the 
volume-preserving systems are so important in statistical physics, but left it 
as a topic for the forthcoming papers. 
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